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The energy flow, dE/dt], is studied at large pseudorapidities in proton-proton colli- 
sions at the LHC, for centre-of-mass energies of 0.9 and 7 TeV. The measurements are 
made in the pseudorapidity range 3.15 < < 4.9, for both minimum-bias events 
and events with at least two high-momentum jets, using the CMS detector. The data 
are compared to various pp Monte Carlo event generators whose theoretical models 
and input parameter values are sensitive to the energy-flow measurements. Inclu- 
sion of multiple-parton interactions in the Monte Carlo event generators is found to 
improve the description of the energy-flow measurements. 
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1 Introduction 

Proton-proton collisions at the Large Hadron Collider (LHC) allow Quantum Chromodynam- 
ics (QCD) processes to be investigated through the measurement of the average energy per 
event (energy flow) in specific angular regions. Such a measurement is useful in examining the 
complex final states that result from a hadron-hadron interaction. Exploiting the large calor- 
imeter coverage of the Compact Muon Solenoid (CMS) detector allows the energy flow to be 
measured over a wider range than was accessible in previous analyses. 

At the LHC, the accessible fraction x of the proton momentum carried by partons can become 
very small. At small x, parton densities become large, leading to an increased probability for 
more than one partonic interaction. The final state resulting from a proton-proton interaction 
can be described as the superposition of several contributions: the hard interaction, initial-state 
and final-state radiation, hadrons produced in additional multiple-parton interactions m, and 
beam-beam remnants resulting from the hadronisation of the partons that did not participate 
in the hard scatter. The underlying event (UE) is defined as everything except the hard inter- 
action, i.e., multiple-parton interactions, beam-beam remnants, and initial-state and final-state 
radiation. The UE plays an especially important role at high centre-of-mass energy. 

Multiple-parton interactions are not well understood theoretically and a systematic descrip- 
tion in QCD remains challenging [2J. Phenomenological approaches to multi-parton dynamics 
rely strongly on parameterised models and tuned parameters (tunes) to describe data. Mea- 
surements of the UE structure have been performed for central values of the pseudorapidity 
[r]), where // = — In [tan (0/2)], with Q being the polar angle of the particles with respect to 
the beam axis 13 -5J. The energy flow has been measured previously in pp collisions [6] at 
= 0.2 — 0.9 TeV, as well as in ep collisions [7J. The extension of the measurements to large 
pseudorapidities and higher centre-of-mass energies is a challenge for the models since, in this 
region of phase space, parton showers (initial-state radiation), as well as multiple-parton inter- 
actions, are expected to play a significant role. 

In this paper, the measurement of the energy flow in the pseudorapidity range 3.15 <\rj\ < 4.9 
is presented. The energy flow in an event is defined as dE/di] = £, /A?/, where is the 
summed energy of all charged and neutral particles in the event measured in bins of pseudora- 
pidity, and A}] is the bin width in r]. Two different centre-of-mass energies, 0.9 and 7 TeV, were 
investigated for two different event classes: minimum-bias events and events with central high- 
transverse-momenta dijets. The latter are expected to be more sensitive to perturbative QCD 
phenomena. 

This paper is organised as follows: Section |2] describes the experimental apparatus, and the 
event selection is explained in Section |3] The main features of the Monte Carlo (MC) event 
generators used for the comparison with data are presented in Section |4] Sections |5] and |6] 
discuss the corrections and systematic uncertainties related to the energy-flow measurement. 
The results are presented in Section|7| and the final conclusions are summarised in Section|8[ 

2 Experimental Apparatus 

This section briefly summarises some features of the CMS apparatus relevant for the present 
measurement. A complete description of the CMS detector can be found elsewhere |8|. The 
CMS experiment uses a right-handed coordinate system, with the origin at the nominal inter- 
action point (IP), the x axis pointing towards the centre of the LHC ring, the y axis pointing 
upwards, and the z axis along the anticlockwise-beam direction. The azimuthal angle (p is mea- 
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2 Experimental Apparatus 



sured with respect to the x axis, in the x — y plane, and the polar angle 9 is defined with respect 
to the z axis. 

The central feature of CMS is a superconducting solenoid, of 6 m internal diameter, providing 
a magnetic field of 3.8 T. Within the field volume are the silicon pixel and strip tracker, the 
crystal electromagnetic calorimeter and the brass/ scintillator hadron calorimeter. The tracker 
measures charged particle trajectories in the pseudorapidity range |?/| < 2.5. 

For triggering purposes, the CMS trigger system ||9l[l0l was used, together with two elements 
of the CMS detector monitoring system, the beam scintillation counters (BSC) IITTl and the 
beam Pick-up-timing for the experiments (BPTX) devices |12|. A BSC detector is located on 
each side of the interaction region, at a distance of 10.86 m from the nominal crossing point, 
covering the \rj \ range from 3.23 to 4.65. Each BSC consists of a set of 16 scintillator tiles. They 
provide information on hits and coincidence signals with an average detection efficiency of 
96.3% for minimum-ionising particles and a time resolution of 3 ns, compared to a minimum 
inter-bunch spacing of 50 ns for collisions. Located around the beam pipe at distances of 175 m 
on either side of the IP, the two BPTX devices are designed to provide precise information on 
the structure and timing of the LHC beams, with a time resolution better than 0.2 ns. 

For the present analysis, the energy was measured with the hadronic forward (HF) calorime- 
ters, which cover the region 2.9< \rj\ <5.2. The front face of each calorimeter is located 11.2 m 
from the interaction point, at each end of CMS. The HF calorimeters consist of iron absorbers 
and embedded radiation-hard quartz fibres, providing a fast collection of the Cherenkov light. 
The collected light is detected using radiation-hard photomultiplier tubes (PMT). The quartz 
fibres are grouped into two categories, segmenting the HF calorimeters longitudinally: long 
fibres, which run over the full depth of the absorber (165 cm lOA/), and short fibres, which 
start at a depth of 22 cm from the front of the detector. Each set of fibres is read out separately. 
The different lengths of the fibre sets make it possible to distinguish between electromagnetic 
and hadronic showers. 

Calorimeter cells are formed by grouping bundles of fibres. Clusters of these cells (e.g., 3x3 
grouping) form a calorimeter tower. There are 13 towers in //, each with a size given by A?/ 
0.175, except for the lowest- and highest-|7/| towers with A/y 0.1 and A?/ ?a 0.3, respectively. 
The azimuthal segmentation of all towers is 10°, except for the one at highest-|^|, which has 
h<p =20°. 

Corrections were applied to the data in order to account for geometrical non uniformities of 
the HF calorimeters, due to non sensitive areas. These are caused by mechanical structures of 
the 20° ^-segment assembly, which are not included in the detector simulation. The relevant 
correction factors for the non uniformities range between 0.98 for the innermost (large tj) pseu- 
dorapidity segment and 0.90 for the outermost (small //) segment, and were applied to the total 
energy deposited in each tower. 

The calibration of the HF calorimeters is based on test-beam measurements with 100 GeV pi- 
ons and electrons [13J. For the relative calibration, a ^"^Co radioactive source was used. The 
energy-flow analysis was restricted to the pseudorapidity range 3.15< \rj\ <4.9. The lowest 
and highest | rj \ towers were excluded from the measurement since they were not properly de- 
scribed in the geometry used for the detector simulation. An additional ring of towers was 
removed from the analysis because it is partially located in the shadow of the electromagnetic 
endcap calorimeter. 
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3 Event Selection 

The data were collected during the first data-taking period of CMS in 2010, corresponding to 
integrated luminosities of 239 }ib^^ and 206 jib^^ for the ^/s = 0.9 and 7 TeV samples, respec- 
tively. In less than 1% of the events more than one interaction vertex was recorded (pile up). 
The energy was measured with the HF calorimeters by summing all energy deposits in the HF 
towers above a threshold of 4 GeV. This threshold, determined from non collision events, was 
chosen to suppress electronic noise. In addition, events in which particles hit the photomul- 
tipliers and cause large signals in the HF calorimeter towers were removed from the analysis 
by means of dedicated algorithms. The rejection criteria are based on the topology of energy 
deposits and the pulse shape/ timing of the signals in the HF calorimeters. 

3.1 Minimum-bias Events 

Minimum-bias events were selected by requiring a trigger signal from the two BPTX detectors 
indicating the presence of both beams crossing the interaction point, coincident with a signal in 
each of the BSC detectors. At least one reconstructed primary vertex was required, with a z co- 
ordinate within 15 cm of the centre of the beam collision region. The vertex was reconstructed 
in a fit with a minimum of four associated tracks [14]. Beam-induced background events pro- 
ducing an anomalously large number of pixel hits were rejected by requiring that at least 25% 
of the tracks in events with more than 10 tracks were of high quality. Tracks were considered 
of high quality when several requirements were fulfilled, such as the normalised X'^ of the fit, 
the compatibility with the primary vertex, and the number of hit layers. A detailed description 
of the high quality criteria can be found in Ref. 1,151 . 

A large number of single-diffraction events was removed with the BSC trigger requirement. 
The fraction of events from single- and double-diffractive dissociation remaining in the samples 
was estimated with a PYTHIA6 MC event sample: 5% (3%) for ^/s = 0.9 TeV (7 TeV). After 
applying all the selection criteria, 2.8 x 10^ and 9.4 x 10^ events remained in the ^/s = 0.9 TeV 
and \/s = 7 TeV data samples, respectively. 

3.2 Dijet Events 

The jets were reconstructed using the anti-A;^ jet clustering algorithm [|T6|, with a distance 
parameter of 0.5. Particle flow objects [fl7| that combine information from all the CMS sub- 
detectors were used as input. The reconstructed jets required only small additional energy cor- 
rections derived from the PYTHIA6 MC generator and in-situ measurements with photon+jet 
and dijet events [18|. The contribution of miss-reconstructed jets was kept at a negligible level 
by means of the procedures discussed in Refs. IIT9ll20l . 

The dijet event sample used in the analysis was a subset of the minimum-bias sample. For the 
event selection, the two highest-pr jets were required to satisfy the condition \A(p{jeti,jet2) — 
n\ < 1.0, where A(p[jet-[,iet2) is the difference in the azimuthal angles of the two jets, and to 
lie within the central region {\ri\ < 2.5). At ^/s = 0.9 TeV (^/s = 7 TeV), the leading and sub- 
leading jets were both required to have pT,jet > 8 GeV (pr,;et > 20 GeV). The requirement on 
rj ensures that the jets are contained in the central region of CMS, outside the acceptance of 
the HF calorimeters. The requirement on the minimum transverse momentum of the jets was 
chosen to reconstruct the jets fully efficiently. The threshold choices result in an approximately 
similar lower limit on the fractional momentum carried by the jets, at both energies. The dijet 
selection criteria resulted in 1.1 x 10^ (1.2 x 10^) events in the ^/s = 0.9 TeV (a/s = 7 TeV) sample. 
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4 Monte Carlo Models 

In this section, the main features of the MC models used for comparison with the data are 
presented, focusing on the implementation and tuning of the UE. Several tunes of the PYTHIA6 
(version 6.420) \\21l event generator were used, each one providing a different description of 
the non diffractive component: D6T [^IJ, DW [22J, Pro-Q20 [23J, Pro-pTO [23J, Z2, as well as 
the central Perugia 2011 tune (Pll) IH and the ATLAS minimum-bias tune 1 (AMBTl) Il25l . 

The tunes differ in the choice of flavour, fragmentation, and UE parameters. The latter set of 
parameters, which are expected to be important for these measurements, includes parameters 
for the parton showers, cut-off values for the multiple-parton interactions, parameters deter- 
mining the geometrical overlap between the incoming protons, and probabilities for colour 
reconnection. Some of the tunes have common flavour and fragmentation parameters, which 
have been determined using data from LEP. With the exception of the parameter settings in 
Z2, Pll, and AMBTl, which are based on measurements from the LHC, the other PYTHIA6 
UE tunes were determined using data from the Tevatron. Tune Z2 is almost identical to tune 
Zl (261, only differing in the choice of parton distribution functions. An overview of the tunes 
can be found in Ref. ||24| . 

The Pll and Z2 tunes, as well as pythiaS ||27|| , use a new model Il28l where multiple-parton 
interactions are interleaved with parton showering. The parton distribution functions used 
are the CTEQ6L [29] set for D6T and Z2, and the CTEQ5L [30J set for the remaining tunes. 
Hadronisation in PYTHIA is based on the Lund string fragmentation model [31J. Predictions 
obtained from PYTHIA8 correspond to the default version of the generator, i.e.,without prior 
tuning to data. In contrast to PYTHIA6, in PYTHIA8 the multiple-parton interactions, initial- 
state and final-state radiation are all interleaved, and hard diffraction has been included. 

The HERWIG++ (version 2.5) [32J MC event generator was also used for comparisons with 
data. It is based on matrix-element calculations similar to those used in pythia. In both gener- 
ators, the evolution of the parton distributions functions with momentum scale is driven by the 
DGLAP ||33fl36l equations. However, HERWIG++ features angular-ordered parton showers and 
uses cluster fragmentation for the hadronisation. HERWIG++ has been tuned independently for 
different centre-of-mass energies. Parameters for the UE and colour reconnections were tuned 
to UE and minimum-bias data at = 0.9 TeV (tune MU900-1) and to UE data at = 7 TeV 
(tune UE7-1) [371. 

Predictions from the event generator CASCADE ||38l |39| are compared to data from the dijet 
event samples. In contrast to PYTHIA6 and HERWIG++, CASCADE is based on the CCFM I|40ti43| 
evolution equation for the initial-state cascade, supplemented with off-shell matrix elements 
for the hard scattering. Multiple-parton interactions are not implemented in CASCADE. 

The DIPSY MC event generator 1441 is based on a dipole picture of BFKL ||45H47l evolution. It 
includes multiple-parton interactions and can be used to predict non diffractive final states. 
Currently, DIPSY is not tuned to experimental data. 

Data are also compared to predictions obtained from pp MC event generators used in cosmic- 
ray physics |48|. EPOSl.99 ||49|, QGSJETlI |50|, QGSJEXOl |51|, and SIBYLL 1521 are considered. 
All models take into consideration contributions from both soft- and hard-parton dynamics. 
In general the soft component is described in terms of the exchange of virtual quasi-particle 
states (a Pomeron at high energies) as in Gribov's Reggeon field theory II53I . At higher energies 
and scales, the interaction is described by perturbative QCD (with DGLAP evolution) by an 
extension of the approach used to describe the soft region. The models were not tuned to LHC 
data. 
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5 Corrections 

The data were corrected to the stable-particle (t > 10^^^ s) level with the help of bin-by-bin 
correction factors derived from simulated events generated with the PYTHIA6.4 MC event gen- 
erator and passed through the CMS detector simulation based on Geant4 f54'|. The position 
and width of the beam spot in the simulation were adjusted to that determined from the data. 
The simulated events were processed using the same analysis chain as for collision data. Cor- 
rections to the measured data were applied bin-by-bin to account for acceptance, inefficiency 
and bin-to-bin migrations due to the detector resolution. 

The correction factor for each bin was determined as the average of the correction factors from 
the PYTHIA6 PO [24J D6T, DW, Z2, and Pro-Q20 tune predictions. All the tunes described rea- 
sonably well the shape of the data before the correction. The deviation from the average was 
included in the model-dependent systematic uncertainty. The number of selected events in 
each of the various MC samples were comparable with the data set size. 

The correction factors were calculated as the ratio of MC predictions at the stable-particle level 
and the detector level. At the stable-particle level, the selected particles were required to be in 
the same pseudorapidity range as that used for the measurements (3.15< \rj\ <4.9), without 
any further requirement on their energies. Neutrinos and muons were excluded. In addition, 
at least one charged particle was required, on each side, within the acceptance of the BSC 
(3.9 < |f/| < 4.4). This requirement was imposed in order to replicate the use of the BSC 
triggers in the detector-level analysis chain. 

For the analysis of dijet events, the kinematic selection of the dijet system was the same for 
the stable-particle level and detector-level jets, i.e., l^ljetl < 2.5, \A(p{jeti,jet2) — 7r| < 1.0, and 
PT,jet > 8 GeV ipT,jet > 20 GeV) for the a/s = 0.9 TeV (^/s = 7 TeV) analysis. 

The event selection criteria applied at the stable-particle level in the MC simulation are sum- 
marized in Table [U 

The bin-by-bin correction factors range from 1.6 to 2.0 (1.5 to 2.0) for the \/s = 0.9 TeV (i/s = 
7 TeV) minimum-bias and dijet data in the four lowet bins of the measurement. Because of 
the larger amount of dead material within the range of the highest- 1 rj \ bin, the correction factor 
is larger, 2.5 and 2.2, for the \/s = 0.9 TeV and ^Js = 1 TeV measurements, respectively. The 
nonlinear behaviour of the HP calorimeters was implemented in the detector simulation and 
was therefore taken into account in the bin-by-bin correction. 

Table 1: Event selection criteria applied at the stable-particle level in the MC simulation. 

Minimum-bias event selection 
N^charged particles > in 3.23 < ?/ < 4.65 and — 4.65 <ri < —3.23 

Dijet event selection 
Ncharged particles > in 3.23 < // < 4.65 and — 4.65 <ri < —3.23 

\n]et\ < 2.5 

\A(p{ieti,iet2) -tt\ < 1.0 
pT,jet > 8 GeV (v/s = 0.9 TeV) 
PT,jet > 20 GeV (y/^ = 7 TeV) 
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6 Systematic Uncertainties 

The systematic uncertainties on the energy-flow measurement are summarised in Table |2] The 
total uncertainty for each \ri\ bin was obtained by adding all the uncertainties in quadrature. 
Depending on the centre-of-mass energy and the event selection, the total systematic uncer- 
tainty was found to be 11 - 15% and 13 - 22% for the minimum-bias and dijet analyses, respec- 
tively. An important systematic effect in the measurement of the forward energy flow is the 
global energy-scale uncertainty of the HF calorimeters, which was estimated to be 10% using 
Z — ?► e+e^ events with one electron in the HF calorimeter and the other reconstructed in other 
subdetectors. 

To estimate the effect of channel-to-channel miscalibration of the HF calorimeters, the response 
per channel was randomly varied between ±10%. The resulting energy flow was shifted by less 
than 1%. To estimate the possible influence of any remaining calorimeter noise, the requirement 
on the minimum energy deposit in a tower was increased to 4.5 GeV. This resulted in a change 
in the energy-flow distributions by less than 2%. 

Variations of the algorithms used to reject events with abnormal signals in the HF calorimeters, 
due to particles hitting the photomultipliers, give a systematic uncertainty on the energy flow 
of approximately 3%. 

The pseudorapidity i] is defined with respect to the origin of the CMS reference coordinate 
system. For events with a primary vertex far (e.g., |z| > 15 cm) from that point, the distri- 
butions of measured variables as a function of ?/ are shifted. To estimate the influence of this 
effect, the energy flow was calculated separately for events with the primary vertex restricted 
to the ranges |z| < 4 cm, 4 < |z| < 9 cm and 9 < |z| < 15 cm. The largest difference was 
approximately 1%. 

The systematic uncertainties due to the model dependence of the bin-by-bin corrections were 
estimated from the differences between the correction factors obtained from the various MC 
tunes and the average value. The differences were applied to the data points uniformly. De- 
pending on the centre-of-mass energy and the | tj \ bin, the resulting changes were between 1 and 
17%. Possible residual systematic effects resulting from the non linearity of the HF calorime- 
ters arise from the difference in the particle spectra and energy distributions predicted by the 
various MC tunes. This was also included in the model dependence of the correction factor. 



Table 2: Systematic uncertainties on the energy-flow measurement for the minimum-bias and 
dijet analyses and their sum. The ranges indicate the variation of the uncertainty within the 
different | rj \ bins. 



Effect 


Minimum-bias analysis 


Dijet analysis 


Energy scale 


10% 


10% 


Channel-to-channel miscalibration 


1% 


1% 


Minimum energy in calorimeter towers 


2% 


2% 


Photomultiplier hits rejection algorithm 


3% 


3% 


Primary vertex z position 


1% 


1% 


Model uncertainty, ^s= 0.9 / 7 TeV 


1-3% / 1-2% 


4-11% / 12-17% 


Non uniformity corrections 


0-3% 


0-3% 


Short-fibre response, = 0.9 / 7 TeV 


3-9% / 3-6% 


6-18% / 6-8% 


Jet energy scale 




2% 


Total, ^/s = 0.9 TeV (7 TeV) 


11-15% (12-13%) 


13-24% (17-22%) 
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The uncertainty on the non uniformity corrections leads to a systematic uncertainty in the en- 
ergy flow of approximately 3% for the highest | rj \ bin, and is negligible for the lowest bin. 

An additional systematic uncertainty comes from the simulation of the short-fibre response in 
the HF, determined by repeating the full analyses using the long fibres exclusively. The former 
was not fully described in the HF simulation at low energies. The resulting uncertainty on the 
energy-flow measurement was between 3 and 9% for the minimum-bias datasets. For the dijet 
samples the uncertainty was between 6 and 10%, except for the two highest-pseudorapidity 
bins for ^/s = 0.9 TeV, where it was between 13 and 18%. 

In the case of the dijet data sample, an additional uncertainty arises from the jet energy scale. 
This was estimated by varying the jet energy by ±10%, leading to an uncertainty in the energy 
flow of 2%. 

The contribution of beam-gas and non interaction events to the event sample was investigated 
by performing the minimum-bias selection on events triggered when there was no beam cross- 
ing. It was found that no such event passed the selection criteria. Other sources of systematic 
effects such as pile up and diffractive modelling were each found to contribute less than 1% to 
the total systematic uncertainty. 



7 Results 

The energy flow was measured with the CMS HF calorimeters at large pseudorapidities, 3.15 < 
< 4.9, and corrected to the stable-particle level. The results are shown in Fig.[l]for minimum- 
bias and dijet events, at ^ys = 0.9 TeV and a/s = 7 TeV. The systematic uncertainties are indi- 
cated as error bars; they are correlated between the bins. The statistical errors are negligible. 
The data are also presented in Tables |3] and |4j 

We observe three distinct features of the data. The first is that the energy flow in both minimum- 
bias and dijet events increases with pseudorapidity, and the increase is found to be steeper for 
minimum-bias events. The second is that the energy flow increases with centre-of-mass energy, 
being a factor of two to three higher at ^/s = 7 TeV than at \/s = 0.9 TeV. The increase in energy 
flow for minimum-bias events is larger than what is observed in the charged-particle multi- 
plicity, dNch/drj, reported in Ref. [55J. Finally, the average energy flow is significantly higher 
in dijet events than in the minimum-bias sample. 

The data are compared to different MC predictions. In Fig. [l] (upper) the minimum-bias mea- 
surements are compared to the predictions of various PYTHIA6 tunes. An interesting observa- 
tion is the similarity between the predictions of tunes Z2 and AMBTl, both based on LHC data, 
and the older D6T tune. Predictions obtained from several other pythia tunes are also shown. 
The variation of the predictions is on the order of 10-20%. 

In Fig. |2] (upper) the minimum-bias measurements are compared to results from different Monte 
Carlo event generators. The PYTHIA6 tunes (also shown in Fig. [l]| are presented as a band, 
which is constructed from the minimum and maximum values of the different PYTHIA6 pre- 
dictions in each bin. The tunes giving the limits of the bands are different for each bin. Also 
shown are the predictions of PYTHIA8, HERWIG++,DIPSY, and PYTHIA6 D6T without multiple- 
parton interactions. The prediction from PYTHIA without multiple-parton interactions is at 
least 40% lower than the measurement. We checked that the same prediction at the parton 
level (without hadronisation) gives smaller energy flow, by approximately 20%. We also found 
that changing the maximum value for the scale used in the parton shower by a factor of four 
increases the energy flow only by ~5%. Therefore, these effects cannot bring the prediction 
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Figure 1: Energy flow as a function of rj for minimum-bias (upper) and dijet (lower) events 
at ^/s = 0.9 and 7 TeV. The data are shown as points with error bars, while the histograms 
correspond to predictions obtained from various PYTHIA6 tunes. The error bars represent the 
systematic uncertainties, which are strongly correlated between the bins. The statistical uncer- 
tainties are negligible. The lower panels show the ratio of MC prediction to data. 
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Figure 2: Energy flow as a functior\ of tj for minimum-bias (upper) and dijet (lower) events 
at i/s = 0.9 TeV and a/s = 7 TeV. The data are shown as points with error bars, while the 
histograms correspond to predictions obtained from various Monte Carlo event generators. 
The yellow bands illustrate the spread of the predictions from the different PYTHIA6 tunes 
considered. The bands are obtained by taking the minimum and maximum variations of the 
PYTHIA6 tunes shown in fig-[l] The predictions from HERWIG++ are made with tunes specific 
to the respective centre-of-mass energy. The error bars represent the systematic uncertainties, 
which are strongly correlated between the bins. The statistical uncertainties are negligible. The 
lower panels show the ratio of MC prediction to data. 
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without multiple-parton interactions into agreement with the data. The HERWIG++ tunes de- 
scribe the measurements well at both centre-of-mass energies. The PYTHIA8 predictions are 
always within the tune uncertainty band of PYTHIA6 and give a slightly flatter energy flow 
distribution than the data. DIPSY, without prior tuning, describes the minimum-bias data well 
for ^Js = 7 TeV. However, it overestimates the energy flow for a/s = 0.9 TeV, by up to 50%. In 
Fig. |3] the measured energy flow is compared to predictions derived from the event generators 
EPOS, QGSJETlI, QGSJETOI, and SIBYLL, which are used in cosmic-ray air shower simulations. 
The description of the minimum bias data by all these models is good. 

The results for the energy flow in dijet events are shown in Figs, [l]-^ (lower). The data are the 
same in the three figures. In general, the rise of the average energy with increasing pseudo- 
rapidity is described by all the models, and the agreement with the data is better than in the 
case of minimum-bias events. The predictions from PYTHIA with different parameter settings 
for multiple-parton interactions agree with the data, whereas the prediction without multiple- 
parton interactions is too low. The prediction from CASCADE, which uses a different parton 
shower evolution, is also too low compared to the data, although slightly higher than the pre- 
diction from PYTHIA without multiple-parton interactions. The energy flow obtained from 
HERWIG++ is in good agreement with the measurement, while that from DIPSY is consistent 
with the data at a/s = 7 TeV, but too high at a/s = 0.9 TeV. The predictions from the cosmic- 
ray MC generators are very close to the data, as shown in Fig. |3] (lower). However, a larger 
deviation is observed in the dijet measurement at a/s = 0.9 TeV, where QGSJETlI and SIBYLL 
underestimate the data in the lowest \rj\ bins, while QGSJETOI is in agreement with the mea- 
surement. 

We have checked that the energy flow in dijet events, after subtracting the energy flow from 
minimum-bias events, is still significantly larger than that predicted from MC models without 
multiple-parton interactions. The disagreement is found to be at least a factor of three. This 
suggests that the energy flow in dijet events is composed of more than a soft underlying event 
and a single parton interaction with a parton shower. 

From the measured energy flow, we can estimate the transverse energy per pseudorapidity 
bin /, given by Ef , = £, sin0,-. The increase of the energy flow with increasing pseudorapidity 
in the minimum-bias events leads to a constant average transverse energy of dEj/ drj ^ 3 GeV 
at ^/s = 0.9 TeV rising to dEj/ drj ^ 6 GeV at ^/s = 7 TeV. The measured transverse energy as 
a function of pseudorapidity is shown in Fig. |4] The error bars were obtained by propagating 
the systematic uncertainties of the energy-flow measurements. 

In dijet events the average energy also increases with increasing pseudorapidity, but the in- 
crease is less steep than in minimum-bias events, leading to a distribution in Ej that decreases 
at a/s = 0.9 TeV from dEj / dy] ^ 4 GeV at |f/| = 3.3 to 2 GeV at |;/| = 4.7. At a/s = 7 TeV 
the transverse energy decreases from dEj/dt] ^ 11.5 GeV at \f]\ = 3.3 to ~ 8 GeV at |?/| = 4.7 
(Fig. |4|. This decrease of transverse energy is consistent with calculations using a pr-ordered 
or virtuality-ordered (Q^-ordered) parton evolution, with the highest transverse parton mo- 
mentum being closest to the hard scatter (at small \rj\), and decreasing towards the direction 
of the proton (large \t]\). The general behaviour of the energy, as well as the transverse energy 
flow, is well described by MC models applying a pj- or virtuality-ordered parton shower with 
multiple-parton interactions. 

The average transverse energy measured in dijet events can also be compared to the transverse 
energy in deep-inelastic scattering (DIS) events measured at HERA [7] at a similar x and Q^. 
We compare our measurement of the transverse energy of dEj / drj ^ 4 GeV at |//| = 3.3 and 
v^s = 0.9 TeV (at a scale of = 4pj ?a 250 GeV-^ and x ~ 0.02) with the corresponding value 
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Figure 3: Energy flow as a function of tj for minimum-bias (upper) and dijet (lower) events 
at ^/s = 0.9 TeV and ^/s = 7 TeV. The data are shown as points with error bars, while the 
histograms correspond to predictions obtained from different cosmic-ray Monte Carlo event 
generators. The error bars represent the systematic uncertainties, which are strongly correlated 
between the bins. The statistical uncertainties are negligible. The lower panels show the ratio 
of MC prediction to data. 
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Figure 4: Transverse energy flow as a function of t] for minimum-bias and dijet events at = 
0.9 TeV and a/s = 7 TeV. The data are shown as points with error bars. The error bars represent 
the systematic uncertainties, which are strongly correlated between the bins. The statistical 
uncertainties are negligible. 



from DIS of dEr/dri ^ 2 GeV at x ~ 0.01 and ~ 250 GeV^. For the measurement at 7 TeV, 
at a scale of = 4p ^ = 1600 GeV-^ and x ~ 0.006, no corresponding measurements from DIS 
exist. This comparison shows that the (transverse) energy flow in pp collisions is significantly 
larger than that measured in DIS, where the contribution from multiple-parton interactions are 
negligible. 

In summary, the shape of the energy-flow distribution, both in minimum-bias and dijet events, 
is reproduced by all the MC event generators that include a contribution from multiple-parton 
interactions. However, the magnitude of the average energy density depends significantly on 
the parameter settings of the different MC tunes. The measured energy flow can thus be used 
for further constraints on the modeling of multiple-parton interactions. 

Table 3: Corrected energy flow dE/drj and systematic uncertainties Ssys for the minimum-bias 
measurements. All values are in GeV. The statistical uncertainties are less than 0.1% in all bins, 
and are therefore not listed. 





Minimum-bias data 
^/s = 0.9 TeV ^/s = 7 TeV 


\v\ 


dE/drj Ssys 


dE/df] Ssys 


3.2-3.5 

3.5- 3.9 
3.9-4.2 
4.2-4.6 

4.6- 4.9 


41 5 
52 6 
78 9 
89 11 
101 14 


90 10 
119 14 
170 20 
220 26 
310 40 



8 Conclusions 

The energy flow at large pseudorapidities, 3.15 < \ri\ < 4.9, has been measured in pp collisions 
for minimum-bias events and events with a dijet system in the central region, \r]jet\ < 2.5, with 
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Table 4: Corrected energy flow dE/drj and systematic uncertainties Sgys for the dijet measure- 
ments. All values are in GeV. The statistical errors are less than 0.1% in all bins, and are there- 
fore not listed. 





Dijet data 
= 0.9 TeV y/s = 7 TeV 


\v\ 


dE/df] 


^sys 


dE/df] Ssys 


3.2 - 3.5 


61 


8 


161 34 


3.5 - 3.9 


73 


10 


206 42 


3.9 - 4.2 


94 


17 


271 49 


4.2-4.6 


101 


19 


348 64 


4.6 - 4.9 


113 


25 


463 85 



PT,jet > 8 GeV {pT,jet > 20 GeV) in pp collisions at ^/s = 0.9 TeV {y/s = 7 TeV). 

By requiring a high-momentum dijet in the central region, the deposited energy increases in 
the entire phase space. Thus the forward energy flow is higher in the dijet data sample than 
in the minimum-bias sample. The increase of the energy flow with increasing centre-of-mass 
energy is reproduced in general by all the Monte Carlo event generators, both for minimiuxi- 
bias and dijet events. The results indicate that predictions of models including multiple-parton 
interactions are close to the data. The Monte Carlo predictions without multiple-parton in- 
teractions, derived from PYTHIA6 and CASCADE, significantly underestimate the energy flow 
compared to data. None of the PYTHIA6 tunes imder study can describe all four energy-flow 
measurements equally well and, in general, they predict a flatter energy-flow distribution in 
minimum-bias events. We observe that the D6T and Pro-Q20 tunes provide the best descrip- 
tion of the minimum-bias and dijet data, respectively. The predictions from the HERWIG++ 
tunes as well as DIPSY are also in agreement with the data. The predictions from cosmic-ray 
interaction models provide the best descriptions of the measured energy flow. 

The variation of the energy flow with rj, both in minimum-bias and dijet events, is reason- 
ably well reproduced by all Monte Carlo event generators with multiple-parton interactions 
included. However, the magnitude of the average energy strongly depends on the parameter 
settings of the different MC tunes, as shown by the large spread in the theoretical predictions. 

A comparison with measurements at HERA in deep-inelastic ep collisions, where the contribu- 
tion from multiple-parton interactions is negligible, shows that the (transverse) energy flow in 
pp collisions is significantly larger. 
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